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A mean field theory of the localization transition for 
bosonic systems is developed. Localization is shown to be 
sensitive to the distribution of the random site energies. It 
occurs in the presence of a triangular distribution, but not a 
uniform one. The inverse participation ratio, the single site 
Green's function, the superfluid order parameter and the cor- 
responding susceptibility are calculated, and the appropriate 
exponents determined. All of these quantities indicate the 
presence of a new phase, which can be identified as the Bose- 
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The localization transition in disordered systems has 
been a focus of statistical and condensed matter physics 
ever since the famous paper of Anderson [lj. Though 
most of the work deals with the metal-insulator transition 
in fermionic systems, there is a recent surge of interest 
in bosonic models pjj^j. The superfluid-insulator transi- 
tion in granular superconductors || or 4 He in disordered 
media J5| are the paradigmatic experimental realizations 
of such systems. The theoretical understanding of phase 
transitions is typically based on a mean field description 
and subsequent fluctuation analysis. The generally held 
belief even today is that the spatial homogeneity of this 
"mean field" allows only for extended states, thus obliter- 
ating the localized phase Q| . Below we demonstrate that 
in a model of hard core bosons with random site ener- 
gies and infinite range hopping this conventional wisdom 
does not hold, and with a suitable choice of the disor- 
der distribution the localization transition can indeed be 
captured within the mean field technique. 

To highlight the involved physics, consider the stan- 
dard argument against a transition in infinite connectiv- 
ity lattices. It is believed that because of the infi- 
nite number of neighbors, every site will be connected to 
"virtually degenerate" sites for a continuous distribution 
of the disordered site energies. Thus hopping between 
these sites always gains kinetic energy with zero cost in 
potential energy, delocalizing the particles. However let 
us observe that for a finite system of size N, the gain 
from hopping between two sites is 0(1/ N) and the po- 
tential energy difference to the energetically closest sites 
is of the same order, since one has N site energies chosen 
independently from a finite interval. This means that 
there is a finite probability that no sites at all are avail- 
able within the 1/N energy window set by the kinetic 
energy. When hopping from such a site, the potential 
energy cost certainly outweighs the kinetic energy gain. 



A complex sum of gains and costs will decide whether a 
state will be localized or extended. 

To demonstrate the nontriviality of this sum, consider 
the spectrum of the ordered array. The ground state 
is homogeneous, non-degenerate with an energy of -1, 
clearly maximally benefitting from the kinetic term. All 
of the other N-l single particle excited states possess zero 
energy. This is so because in this model to ensure or- 
thogonality to the ground state the wavefunctions of the 
excited states have fluctuating signs across the sample. 
This destructive quantum interference frustrates the ki- 
netic term, preventing any gain from the hopping pro- 
cess. So for weak disorder there is more kinetic energy 
to be gained by staying extended, and only a sufficiently 
strong disorder can localize the ground state. For the 
excited states however staying localized at a suitable site 
offers the lowering of the total energy in the absence of a 
kinetic energy premium, thus the excited states will be- 
come localized for arbitrarily small disorder. As we will 
see, this physics can be brought out by a non-traditional 
choice of the disorder: we will focus on the triangular 
distribution of the site energies. 

To establish the framework of the physics, consider the 
Hamiltonian 
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where a\ (a^) creates (annihilates) a hard-core boson at 
site i (i = 1, . . . , N). We consider the case of infinite 
range hopping: Jy = iV _1 for all pairs, where the mean- 
field approach is exact. J a = /i + hi, where \i is the 
chemical potential and hi is a random on-site energy. 
We recall that Q in the absence of disorder there are 
two phases of the model: for generic fillings the bosons 
can propagate. At zero temperature they form a super- 
fluid. On the other hand for precisely one boson per site, 
the particles localize in a Mott insulating phase. When 
we take away a single particle from this insulator, the re- 
sulting hole will behave exactly as the first particle added 
to the empty lattice. This can be proven by performing 
a particle-hole transformation on the Hamiltonian. 

In what follows we approach the localized phase from 
two directions. First, from the Mott insulator by adding 
holes. In this case a one particle approach clearly suf- 
fices. We calculate the density of states (DOS) and the 
participation ratio, and show that above a critical disor- 
der the ground state becomes localized. Critical expo- 
nents will be evaluated as well. Second, we approach the 
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same phase by decreasing the density of holes from the 
extended phase and compute the superfluid order param- 
eter and the susceptibility. 

In studying the one particle problem, the quenched dis- 
order averaging < • • • > aV e is performed with the replica 
trick, allowing us to rewrite the DOS as: 



p(A) - <^6(X + J fe ) > ave 

k 
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= lim »-°]v^ Im aA^ J dhp ( h ) 
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A is an energy eigenvalue with a positive infinitesimal 
imaginary part and the J^'s are the eigenvalues of the J 
matrix. Next the off-diagonal terms are decoupled using 
an auxiliary field z a , and the x integrals are performed 
to transform the exponent into: 



S, ii = -~ z l + N ln < ex P [ - \ ^ ( X - h )} ( 3 ) 
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where the averaging over the disordered site energies 
J dhP(h)A(h) is denoted by < A(h) >. In the expec- 
tation value we use the inverse of the number of lattice 
sites TV as a small parameter to carry out an exact ex- 
pansion. We also keep only the terms linear in the replica 
number n to arrive at the exponent: 



S eff = - 



< 



X 



zl + Nn< ln(A -/*)>) 



Finally performing the integration over the auxiliary field 
z a yields: 



P(A) + 1<5(A + A ), 



(4) 



where Ao is obtained from: < (Ao — ^) 1 >= 1- We 
assume that P(h), the distribution of the site energies is 
non-zero on the interval (—A, A). 

The DOS was easiest to obtain in the above path- 
integral framework, but the subsequent physical quan- 
tites can be determined by the simpler method of calcu- 
lating the eigenvectors <p(X) of the J matrix. One finds 
tpi(X) = m/(X — hi), where m is the superfluid order 
parameter: m — (l/7V)Ev«- The self-consistency equa- 
tion for m yields: 1 = 1/iVE 1/(A — hi)- If A is inside 
P{h), then a few ipi(X)'s will be ~ 0(1) and the rest 
of the ipi(X)'s will be ~ 0(1/ N) to satisfy the normal- 
ization condition, i.e. the states inside the continuum 



are localized, whereas for Ao outside the continuum all 
¥>i(Ao) ~ 0{\/\/N), describing an extended state. This 
picture is in complete accordance with the above deter- 
mined DOS. The ground state becomes localized if — Ao 
reaches the bottom of the band. This does not happen 
for a rectangular distribution since explicit calculation 
yields Ao = AcothA > A. However for the triangu- 
lar distribution P(h) = A~ 2 (A — \h\) the ground state 
eigenvalue — Aq is given implicitly by the equation 



A 2 = (A - A) In 



An- A 



+ (A + A) In 



A + A 



(5) 



Ao Ao 
which is valid only for Ao > A. Taking the Ao — > A + 
limit in Eq. (||) shows that — Ao reaches the bottom of 
the band at the critical disorder A c = 2 In 2 = 1.38. For 
A > A c the DOS becomes identical to the distribution 
of the site energies, i.e. p(A) = P(X). 

One measure of localization is the participation ratio 
V = (£, N 2 )7(Ei l^il 4 )' V is proportional to the 



system size N for extended states and remains 
for localized ones. The expression for tfi yields 

V_ = <{X -h) 
N 



-2 > 2 
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(6) 



< (Ao - h) 

for A < A c and V/N = for A > A c . For the uniform 
distribution there is no critical disorder, and the partic- 
ipation ratio is V/N = 3/(3 cosh 2 A + sinh 2 A) ~ 0(1), 
clearly indicating that the ground state remains extended 
for any finite disorder strength. For the triangular dis- 
tribution the result of the integral (||) is 

A2\2 X 2 a2 

* ln 2 ^A, (7) 
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and Ao is obtained from Eq. (|^) for A < A c . As Ao 
approaches A, V/N disappears as oc (A c — A) 2 , indicat- 
ing that for A > A c the ground state becomes localized. 
From now on we will focus only on this more interesting 
case of the triangular distribution. 




FIG. 1. Analytical and numerical results for the ground 
state participation ratio. 
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One can independently test these results numerically 
by diagonalizing the matrix J and measuring V . Fig. 1 
displays convincing agreement between analytic results 
(solid line) and numerical results (symbols). The cross- 
ing and subsequent decrease of V with N (inset) pro- 
vides compelling evidence for the localization of the 
ground state. To demonstrate how V distingishes be- 
tween extended and localized states, one usually con- 
siders two possibilities, ipi = 1/yN and (fii = Si^ . 
These yield V = N and V = 1 respectively. One of- 
ten equivalently phrases this analysis in terms of the in- 
verse participation ratio V^ 1 which takes on small val- 
ues (D(l/N) for extended states and values near unity 
for localized ones. However, if ipi has a few large am- 
plitudes on selected sites, and an extended background, 
tpi = [l/2]5 Mo + l/y/2(N - 1)(1 - 6i, io ), V and V^ 1 will 
be close to one and indicate, incorrectly, that the state 
is localized. To eliminate the possibility of misidcntify- 
ing such "pseudo-localized" states, we computed 'P _1 (n), 
which we define by systematically removing the n sites 
where the wavefunction assumes its largest values. We 
show the results of this calculation in Fig. 2. 
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FIG. 2. The truncated participation ratio. Here the 
four indices refer to the ground state, first excited state, a 
mid-band and the maximal energy wavefunctions. 

When A = 1.0 < A c , the ground state is extended, 
and 'P~ 1 (n) remains 0(1/N) as sites are removed. Mean- 
while, for all the excited states, V~ Y rapidly plunges to 
much smaller values than 0(1/N), emphasizing that as 
a few sites are removed from the sum, the weight of the 
remaining sites is negligible. We checked this behavior on 
different lattice sizes, and found that , P~ 1 (n) fell slightly 
more rapidly with n for larger lattices, eliminating the 
possibility that the states were extended but over a small 
fraction of the lattice. We conclude that all states in the 
continuum are indeed localized. When A = 1.7 > A c 
even the ground state is localized, as shown. 

Returning to the analytic calculations,the transition 



can be parametrized by the chemical potential as well. If 
/i exceeds a critical value /x c , then the ground state has 
no holes. At fi c the first hole appears in the system. [i c 
is given by the lowest eigenvalue of the — J matrix, i.e. 
/i c = Ao for A < A c and \i c = A for A > A c . The 
transition is well captured by the imaginary time on-site 
Green's function B 



7(r) 



i^(r rfl J(r)a i (0)) o = J 



d\p(X) exp[— r(/i — A)], 
(8) 



for t > 0, where (■ ■ -)q is the ground state expecta- 
tion value, and T T is the time ordering operator. For 
weak disorder the transition happens to the state at — Ao. 
Combining Eqs. (Q) and (^|) one obtains at criticality: 
g(r) w X/N + exp[— r(Ao — A)] near the transition, for 
t ^> 1. On the other hand, for strong disorder (A > A c ) 
the critical g(r) decays as: g(r) cx r~ 2 . This change of 
the critical behavior suggests that we enter into different 
phases for A < A c and A > A c . It is worth noting that 
if the distribution P(h) cx (A — h) a , then in the strong 
disorder regime g(r) oc r _a_1 right at the critical point. 
This means that in the strong disorder regime the critical 
exponent depends on the distribution of the disorder. 

If one approaches the same localized phase from higher 
densities of holes, these one particle techniques must 
be abandoned. We will concentrate on the free en- 
ergy of the Hamiltonian ([[]): / = ksTlnZ/N, where 
Z = Tr exp —f3H. Introducing the magnetization m as a 
Hubbard-Stratonovich field decouples the different sites 
to yield a single site problem: 



f = —m 2 lnQ,- 



where 



Qi = Tr exp 0[(p + hi) a^a + ma) + m*a] 



(9) 



(10) 



For hard-core bosons the Hilbert space is only two di- 
mensional, as a site may be only empty or occupied by 
one particle. This allows the exact evaluation of Q and 
hence the free energy per site, which we give here only 
for T = 0: 



1 



1 



f = -m 2 --<(p- h) + v /( A1 -/ l )2 + m 2 > . (ii) 
m is determined by the saddle point condition: m = 



m < [(/i - h) 2 + 



-1/2 ^ xhe order parameter m 
is proportional to (aj). The hard-core boson problem is 
equivalent to a spin 1/2 XY model in a transverse field, 
and m corresponds to the magnetization in the XY plane. 

m = is always a solution, but upon exiting the Mott 
insulator below a critical chemical potential fi c another 
solution with nonzero magnetization m appears. For 
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weak disorder A < A c /i c = Ao- Expanding in m near 
the critical point the magnetization behaves the same 
way as it does without disorder, i.e. m oc y/JT c — /i, the 
well-known Landau result. 

In the strong-disorder regime when increasing the den- 
sity, the particles occupy localized states. Do we have a 
superfluid in this case? In order to get the answer, we 
write out the saddle-point equation for the triangular dis- 
tribution: 



/•*+ m 



\i- ln 



H+ + \Ui% + 



(12) 



H- + \ I H 2 - + rn 2 



H + yV + 



+ 2 a/ (i 2 + m 2 - \J [i\ + rn 2 - \J fJ/L + rn' 1 



where we have introduced the notation fi+ = \i + A and 
//_ = ji — A. For strong disorder fi c = A, because if jj, 
is slightly below A, i.e. /i_ is negative, one sees immedi- 
ately from the second term of the right hand side (RHS) 
of Eq. (13) that m must be different from zero, otherwise 
the logarithm 'blows up'. Clearly a superfluid is formed, 
even though the constituent particles possess localized 
wavefunctions. The physics is then that the phases of 
the localized wavefunctions lock up to yield the off diag- 
onal long range order. This situation is very reminiscent 
to the "localized superconductor" picture of Ma and Lee 



The above results clearly show that the truly localized 
phase occupies only a line in the /i — A plane. However 
we think that this is indeed the seed of the "Bose-glass" 
phase, because for the finite range hopping model if one 
partitions the system into blocks of the size of the hop- 
ping length, each will support only one of these truly 
localized, non-superfluid states. There will be a macro- 
scopic number of these blocks, thus the original line will 
expand into a finite region as a function of the density. 

The one-particle states being localized, it is natural to 
assume that m <C (/i c — //) = |/x_|. In this case the second 
term of the RHS of Eq. (13) behaves like //_ ln(m 2 /|/i_ |), 
and this term must be finite, so we find 



m oc y/ He — H exp ( — 



Hc~ H 



(13) 



where a = A(A — A c )/2. As we can see, below \Xc the 
system is superfluid, though the magnetization is much 
smaller than in the weak-disorder case. The critical be- 
haviour is again different in the weak and strong disorder 
regime. It is remarkable that in the study of the corre- 
sponding one dimensional problem also essential singu- 
larities were found H. 

Finally we calculate the susceptibility by adding an 
infinitesimal in-plane field B to m and taking \ = 
— (d 2 f/dB 2 ). Approaching the transition from the m > 



side for arbitrary disorder we get x oc (/i — /i c ) _1 . On 
the other hand approaching from the non-superfluid side 
(i.e. m — 0), while for weak disorder one again obtains 
X oc (/U — /U c ) , for strong disorder \ remains finite even 
at the transition. This difference of the exponents again 
demonstrates that the transitions from the Mott insula- 
tor into the superfluid or into the localized region belong 
to different universality classes, further strengthening the 
argument that the localized region in fact is a separate 
phase. Just as for the Green's function, the above ex- 
ponents also depend on the asymptotics of the disorder 
distribution for strong disorder. 

To summarize, we investigated the disordered hard 
core boson problem. We proved that, contrary to pre- 
vious beliefs, the mean field theory is able to capture 
the localization transition. We approached the localized 
region both from the Mott insulator and from the super- 
fluid phase, calculating the density of states, the inverse 
participation ratio, the on-site Green's function, the mag- 
netization and the susceptibilty. We tested our theory 
with independent numerical investigations and found de- 
tailed agreement. The important observation is that the 
exponents of all of the above physical quantities differ for 
the direct insulator - superfluid and the insulator - glass 
or superfluid - glass transitions, clearly demonstrating 
that the glass transition belongs to a new universality 
class. We also gave a real space blocking argument why 
we expect the glass to expand into a region in the pa- 
rameter space when the range of hopping is reduced to 
finite values. 
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